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1. Introduction 



In [1], V. Korepin and the author brought up the issue of the sensitivity of the six- 
vertex model to its boundary conditions (even in the thermodynamic limit). The moti- 
vation came mostly from some recent work on domino tilings [2,3,4], in which boundary 
conditions seemed to affect greatly the typical arrangement of dominos. The problem of 
counting domino tilings is equivalent to the six-vertex model with particular Boltzmann 
weights; this is schematically described on Fig. 1. Therefore it seems natural to investigate 
the corresponding problem for the general six-vertex model with arbitrary weights. 




Fig. 1: Correspondence between vertices of the six- vertex model and small 
patches of a domino tiling. 

The usual studies of the six-vertex model (see [5] and references therein) are made 
by assuming periodic boundary conditions (PBC). In [1], different boundary conditions, 
the so-called domain wall boundary conditions (DWBC), were used (Fig. 2a), and the 
thermodynamic limit of the model was investigated using determinant formulae for the 
partition function [6,7]. The main result found was an expression for the bulk free energy 
in the disordered phase of the model, which is different from the usual expression for the 
case of periodic boundary conditions. It should be noted that the DWBC correspond to 
the Aztec shape in the domino tiling language (see Fig. 2), which is precisely the type of 
tiling which was considered in [2,3]. 
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a) 



b) 



Fig. 2: a) A configuration of the six-vertex model with DWBC, and b) one 
possible corresponding tiling of the Aztec diamond. 

Here, we use a new method to compute the bulk free energy with DWBC in all phases 
of the model; in particular, we obtain an independent confirmation of the results of [1]. 
In section 2, starting from the determinant formula for the partition function, we shall 
rewrite the latter as a matrix integral, but with a measure on the space of hermitean 
matrices which is not necessarily smooth. In the disordered phase (section 4), the measure 
will turn out to be smooth, whereas in the ferroelectric and anti-ferroelectric phase (section 
3 and 5) it will be discrete (when expressed in terms of the eigenvalues). The size N of the 
matrices is the size of the original square lattice, and therefore the thermodynamic limit 
can be investigated using tools from large N matrix models. Since the results of section 
5 (concerning the anti-ferroelectric phase) are new, they are analyzed in more detail by 
considering various limits of the parameters, and the subleading correction of the free 
energy is calculated. 

2. Properties of the determinant formula 

We use the same notations as in [1]. We consider the homogeneous six- vertex model, 
with the following parameterization of the Boltzmann weights attached to the vertices: 



a = sinh(t — 7) 



b = sinh(t + 7) 



c = sinh(27) 



(2.1) 



The domain wall boundary conditions (DWBC) mean that external horizontal arrows 
are outgoing, whereas external vertical arrows are incoming (Fig. 2a). These boundary 
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conditions only exist for square lattices. In [6,7], it was shown that the partition function 
of the six-vertex model with DWBC on a N x N lattice could be written as: 

(sinh(t + 7 ) sinh(t- 7 ))" a 
= ; ~2 T N 

N-l 1 ^ 



where tn is a Hankel determinant: 

" ^i+k-2 



tn = det 

l<i,k<N 



fai+k-2 



(2.3) 



Here, 



6(t) = sinh(2 7 ) (2 4) 

[ ' sinh(t + 7) sinh(t - 7) 1 ' ' 



It is known that such determinants are tau-functions of the Toda semi-infinite chain hi- 
erarchy in terms of appropriate parameters. Here, as a function of t, the tn satisfy the 
ususal Toda equations under the bilinear form [8,1]: 



TNT N ~ T' N 2 = TN+lTN-l ViV > 1 (2.5) 

This equation was used in [1] to derive the bulk free energy of the model in the ferroelec- 
tric and disordered phase by making an appropriate Ansatz on the large iV form of tn- 
Unfortunately, the Ansatz in the anti-ferroelectric phase is not that simple, as we shall see, 
and would be hard to justify at this point. 

We shall therefore use another approach here, based on the equivalence of Hankel 
determinants with one-matrix models [9,10]. Let us write formally <f>(t) as a Laplace 
transform: 

<f>(t) = J dm(X)e tx (2.6) 
where dra(A) is a measure. We then notice that the derivatives of <f>(t) are the moments: 

^-<j>(t) = J dm(\)Xe tx (2.7) 

Inserting this into (2.3) leads to: 

r N 

tn = J dm(\ 1 )...dm(\ N ) (-O'll [e Ul A: +fT(l) - 2 ] (2.8) 



We see that appears naturally the Van der Monde determinant A(Aj) = det(A^ 1 ) = 
Yli<j(^i ~ ^j)- After a few elementary manipulations we find: 

™ = W.J dm(Al) • • - dm ( X N)HK) 2 e^ Xi (2.9) 

If dm(A) is a smooth positive measure of the form dm(A) = dAe _1/ ( A ), then we recognize 
in (2.9) the expression in terms of its eigenvalues of the matrix integral: 

t n ~ JdMe tT \- tM + V (M)} (2 . 10) 

where M is a hermitean A x A matrix, and dM is the flat measure. 

As we shall see, if the measure is not smooth, we shall end up with expressions 
which can still be treated using appropriately adapted matrix model techniques. This 
is typically the case of discrete measures that appear in sums over Young diagrams 
[11,12,13,14,15,16,17,18]. 

Expressions of the type (2.9) have been widely studied in the literature (on random 
matrices in particular). One important goal is to find their large iV asymptotic behavior. 
Here we shall mention the simplest method to find their leading large A behavior: the 
saddle point method. The basic idea is that log A(A^) 2 , being a sum of ~ A 2 terms, scales 
as A 2 in the large limit, whereas there are only A variables of integration. Therefore the 
integral is dominated by a saddle point. An important remark is that, in order to find the 
saddle point, we must write our action (i.e. log of the function integrated) in such a way 
that all terms are of the same order A 2 . Here, the term t K is naively of order A, and 
we reach the important conclusion that the Xi will scale as 

\ z oc N(n (2.11) 

After the change of variables Aj — > fa, one can use the saddle point approximation, which 
gives us access to the function / defined by 

/ = hm lQg(TjV / Cjv) (2.12) 

where cn = (Il^=o n ') 2 - / i s essentially the bulk free energy, cf Eq. (2.2). Note that the 
saddle point is a very crude approximation in the sense that it does not naturally allow for 
a systematic computation of subleading corrections; however it will be sufficient for our 
purposes. We now proceed with a separate discussion of the different phases of the model. 
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3. Ferroelectric phase 

This is the phase in which the weights are given by (2.1) with t and 7 real, I7I < t. 
We use the following decomponsition: 

= . _ sinh ( 2 7) = 4 y e -m sinh(2 z) (31) 

rv ; sinh(t + 7) sinh(t - 7) ^ v ' ; v ; 

We are in the situation where the measure dm is discrete. The determinant takes the form 

T/v = 2 N " A (^) 2 e" 2t ^ h I] smh(2 7 /,) (3.2) 

h,...,l N —Q i 

(we have neglected here, as in all subsequent calculations, constant factors which manifestly 
do not contribute to the bulk free energy). This expression is very close to what one 
encounters when studying the Plancherel measure (or other similar measures) on Young 
diagrams [11]. In the context of Young diagrams, the k represent the shifted highest 
weights li = rrii + N — z, where the m; are the usual highest weights (sizes of the rows 
of the diagram), and one is usually interested in the limiting shape of the Young diagram 
when its size is sent to infinity. There has been a lot of work on this type of expressions, 
both in the mathematical literature [11,15,16,17] (the recent work being concerned with 
fluctuations around the limiting shape, which we shall not discuss here) and the physical 
literature [12,13,14,18]. One relevant observation from [12] is the following: after the 
rescaling p = l/N, all sums look like Rieman sums and one is tempted to replace them 
with integrals, and then apply the saddle point method. This is correct on condition that 
one imposes an additional constraint coming from the discreteness of the li. In Eq. (3.2), 
all li must be distinct integers (due to the Van der Monde determinant), and therefore 

\h-lj\>l Vi^j (3.3) 

If we introduce the density p(y) dfi of the \ii = h/N, normalized so that f p(n) dfi = 1, 
then (3.3) implies that it must satisfy the inequality 

P {H) < 1 (3.4) 

In general, when the k are trapped in a well of the potential (as is the case here), there 
will be a saturated region at the bottom of the well where p(p) = 1, and an unsaturated 
region where p((i) < 1. 
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Let us now proceed with the solution. Once the rescaling ^ = l t /N is performed, 
one notices that up to corrections exponentionally small in N, smh(2'jN/ji) xs ^e 2 ' 7 '^ 1 . 
Therefore 

t n *c' n 2 n2 Afafe-W-WEt* (3.5) 

/j, 1 ,...,lj, N €-p r Z + 

where c' N = N N . Of course, once this simplification is made, we regognize a well-known 
expression; in fact, going back now to the original variables ^ one can compute tn directly 
using the Cauchy identity for Schur functions. However, to emphasize the similarity with 
the other phases (which do not possess such a simple group-theoretic interpretation), we 
shall use the saddle point method, following the solution of [18]. Since tn only depends 
on t — |7|, we temporarily set 7 = 0. 

The support of the saddle point density p(/i) is expected to be of the form [0, /?]; the 
saturated region is [0, a], whereas the unsaturated region is [a, 0\. We define the resolvent 

U ^ = / ~7 — V (3 - 6) 

Jo z — [X 

for all complex z £ [0,(3]. The saddle point equations can be written in terms of u: 

u(/j + i0) +u(/j-i0) = 2t Vne[a,0\ (3.7) 

In order to solve the equation, we first remove the logarithmic cut of ui with the redefinition: 
ui{z) = w(z) — log -j^zr^- is analytic everywhere except on [a, (3} and satisfies 

£>(// + iO) + u/(//-i0) = 2t-21og— — (3.8) 

fi — a 

This completely determines it to be: 



r/3-iO j / / 

Co{z)=t-^{z-a){z-(3) / — ^ ^T~Z ( 3 - 9 ) 

f io 2z7r(z - z') y/(z' - a)(z' - (3) z' - a 



After some calculations, we find that 

y/0(z -a) + ^a{z - (3) 



u>(z) = t — 2 log 



y/z(0 - a) 



(3.10) 



The endpoints a and (3 are determined by imposing w(z) ~ j as z — > 00. This gives rise 
to two equations: 

" io S (3.11) 
v^= 1 
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whose solution is: 



, * n , t 

a = coth - p = tanh - 

2 2 



In order to conclude, one expands further the function u(z): 



, , 1 a + p 1 1 

z 4 z 2 



and uses the fact that 



9f n , s a + 3 

i = - 2 <"> = -nr = cotht 



Integrating once and restoring 7, we have the final result 

<=>j — 



sinh(t — I7I) 
which coincides with what was found in [1]. 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



4. Disordered phase 

In this phase, one usually rewrites the weights 



a = sin(7 — t) b = sin(7 + t) c = sin(27) 



(4.1) 



with redefined parameters t and 7, \t\ < 7, and the function <f>(t) = sin(27)/(sin(t 
7) sin(t + 7)); the partition function is then given by 



(sin( 7 + t) sin^-t))^ 2 
= ; —9 t~n 



(4.2) 



with tn still given by (2.3). The Laplace transform is: 

4>(t) 



sin (2-,) f + °° d Xe tX ^fa-W 



sin(7 + t) sin(7 — t) 



sinh ^7r 



(4.3) 



This time the measure is smooth and tn is a matrix integral in the usual sense. 
We must now rescale the variables A*. We choose to define \ii = ^Xi/N . Then: 



/ -N z 
TN = C N 7 



N 



d//i . . .d^ArA(^) 2 J~J 



sinhiV^(^-l) 
sinh N/j,i^ 



(4.4) 



One then simplifies the potential by using: ^ ~ e N ^\. Therefore, 



tjv ~ c' N 7" 



/+oo 
d» 1 ...d m A(» i ) 2 e N ^> i -^ (4. 5) 
-oo 



Note that the matrix integral only depends on the ratio Q = t/^. 

The matrix model (4.5) is fairly simple and can be solved easily in the large N limit 
via the saddle point method. One introduces again the saddle point density of eigenvalues 
p(p) dp, normalized so that J p{p) dp = 1. The support of p(p) is assumed to be a single 
interval [a, 0\ (a < < /?), due to the shape of the potential (single well centered around 
0). The resolvent is defined as before. The saddle point equations read: 

u(p + i0) + u(p-i0) = -C + signGu) VfjLe[a,0\ (4.6) 

where the right hand side is simply the derivative of the potential. The solution of this 
equation: 

a/ f3(z — a) — i\J — ol{z — (3) 



<jj(z) = — — + — log 

2 ITT 



(4.7) 



y/z(J3 - a) 

is very similar to the ferroelectric phase; and the rest of the calculation goes along the 
same lines. 

Requiring that u)(z) ~ \ as z — > oo, we obtain the 2 equations: 



1 _ t - ± loe v^±!^ 

1 ^ - l7 r iQ g (4.8) 

v 7 — a/9 = 7T 



which we solve for a and /?: 



Noting that 



we find 



a = -7rtan^(l - C) /? = tt tan ^(1 + () (4.9) 

f = -logcos^C + cst (4.11) 

We shall not discuss how to fix the constant of integration, since this will be addressed in 

the next section in a more general setting. Reintroducing the 7 dependence coming from 
Eq. (4.5), we have the final expression: 

ef = — -J-p (4.12) 

2 7 cos 1 } 

which reproduces the result of [1] . 
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5. Anti-ferroelectric phase 

We finally study the most interesting phase, in which the weights are given by 



a = sinh(7 — t) b = sinh(7 + 1) c = sinh(27) (5-1) 

with \t\ < 7, and the partition function by 

(sinh( 7 + t) sinh^-t))^ 2 

Z N = -2 TN [0.2) 

' N-l x " 



with <p(t) = sinh(27)/(sinh(7 + t) sinh(7 - £)). 

5.1. Bulk free energy 

We have the expansion 

J>(t) = - Sinh(27 ^ , = 2 V e™e-^ (5.3) 

^ v ; sinh( 7 + 1) sinh( 7 - 1) v ; 

We perform the rescaling \ii = 2^U/N and find that tn takes the form: 

r N = c' Nl - N2 £ A(^e N ^M^^ (5.4) 

jUi,...,MJv€^Z 

The remarkable feature is that Eq. (5.4) is identical to Eq. (4.5) up to the discrete nature 
of the variables! We shall comment on this later. 

The situation is a bit more complicated than in the previous cases, since we now expect 
a saturated region [a', 13'} at the bottom of the well (a' < < (3') and two unsaturated 
regions [a, a'] and [/?',/?] on each side. This is a two-cut situation, which is in fact the 
reason why the naive approach of [1] fails in the anti-ferroelectric phase (see section 5.3 
for more on this). Let us define as before £ = i/7, the density p(pt) and its resolvent u>(/j,). 
The constraint coming from the discreteness of the \ii reads 

P (jjl) < ±- (5.5) 
27 

Therefore we have in the saturated region the equation 

M = ^ (UJ{P " Z0) " ^ + Z0)) = 2^ V ^ G P ' ] (5 - 6) 
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whereas in the unsaturated regions, the saddle point equations are 



wQi + iQ) + w(//-iO) = -C + sign(» V/i G [a, a'] U [j3' , (3] (5.7) 

with C = t/'j. 

We could proceed as in the previous sections; this would lead to a representation 
of w(z) in terms of elliptic integrals. However, this would be fairly cumbersome and we 
proceed instead as follows. Introduce an elliptic parameterization 

^^^-aK^-aof-^ d " = (5.8) 

2 J p y/(z - a)(z - a')(z - ()')(z - (3) 

which corresponds to setting: = sn 2 (u,k) with k = yj (^/"I^)^-"') ' With an 

appropriate choice of path of integration, this maps the complex plane (resp. upper 
half-plane, lower half-plane) onto the rectangle [0, K] x [—iK',iK'\ (resp. [0,K] x [0, iK'], 
[0,K] x [— zK',0]), where K and are the usual complete elliptic integrals of the first 
kind. Similarly, the second sheet of the double covering is mapped onto the other half of the 
torus, which can be chosen to be [-K, 0] x [—iK', iK']. The point of this parameterization 
is that the resolvent u> is now a well-defined function of u. In fact we have the following 
properties: 

(i) The function ui{u) can be extended to a holomorphic function in the whole u plane. 

(ii) The function uj{u) satisfies the following functional relations (for all complex u): 

u(u + 2iK') =u(u) (5.9a) 

7 

u(u + 2K) = u(u) - 2 (5.96) 

<jj{u) + u(-u) = 1 - C (5.9c) 

Eq. (5.9a) is the analytic continuation of Eq. (5.6). Similarly, by combining the analytic 
continuations of the two equations contained in (5.7), one obtains Eqs. (5.96, c). 

(iii) The function u(u) has the following expansion near Uoo = u(z = oo): 

2 

U = 1, a, W% =7r( U - U oc) +O(u- Uoo ) 2 (5.10) 

y/{P' - a){p - a') 
This is a rewriting of the condition uj(z) ~ \ at infinity. 
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Using properties (i) and (ii) (Eqs. (5.9a, 6)), we conclude that £^uj(u) is a doubly 
periodic holomorphic function, and so is a constant. In order to restore the coefficients of 
u>(u) we can use properties (ii) or (iii). We find that 

w(u) = — ^(u-tioo) (5.11) 
A 

plus several conditions relating the different parameters of the problem: 



— - — 

~K ~ 2^ 

y/{P' -a)(f3-a') = 2K (5.126) 

Uoo _ 1 - C 
K ~ 2 



(5.12a) 



(5.12c) 



Relation (5.12a) is particularly interesting since it shows that the elliptic nome q = 
e -irK'/K _ e -7r 2 /2 7 (jgpgn^g on iy on ry ( anc i no t on Q_ Also, the dual nome (under modular 

transformation) q = e -27 is up to a sign the quantum group deformation parameter of the 
model. 

We can rewrite the three conditions in terms of the endpoints; we find 

P-a = 2K- ^ 



sn«oo cnttoo 



P-a' = 2K CUU T (5.13) 
sn-Uoodn-Uoo 

P-fl = 2K CnU °° dnU °° 
sn-Uoo 

In order to completely fix the four endpoints a, a', /?', /?, we need one extra relation; this 
is the equality of chemical potentials in the two unsaturated regions. This relation takes 
the form 

rf>' 

/ (u(/j, + iO)+u(n-iO))dn=(l-C)^ + (1 + C)«' (5-14) 

J a' 

Using the expression (5.11) of co(u), we can rewrite it as 

P'-(P- (3') Sn ^°° Z( Uoo ) = (5.15) 
where Z is Jacobi's Zeta function; this fixes b' to be 

P' = 2KZ( Uoo ) (5.16) 
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The endpoints are now determined by (5.13) and (5.16), supplemented by the value (5.12c) 

Of Woo. 

At this point, we are ready to calculate the free energy. We first rewrite explicitly the 
resolvent (Eq. (5.11)) under the form 

r°° dz' 

lo(z)= / (5.17) 

Jz y/(z' -a){z' -a')(z'-/3')(z'-l3) 

Next we expand it to order 1/z 2 to find 

df a + a > + ( 3> + (3 

d( = 4 (5 - 18) 

which generalizes Eq. (4.10); using some known identities satisfied by Zeta and theta 
functions, we obtain 

df 7T^(7TC/2) 



d( 2 e 2 (ir(/2) 

where we recall that 6*2(2) is 



(5.19) 



e 2 {z) = 2 J2 q (n+1/2)2 cos(2n + l)z (5.20) 

n=0 

There are a variety of ways to find the integration constant. One is to calculate explicitly 
/ (for a particular value of £, e.g. ( = 0) using this matrix model solution, and then 
restore the 7 dependence coming from (5.4); this is a straightforward but tedious exercise. 
Another possibility is to use the known limits £ — > ±1, that is t — > ±7, where we should 
have (see [1]) 

e f — (5.21) 



Either way, we finally find: 



e "2 7 « 2 (|l) (6 ' 22) 



where we recall that the elliptic nome is q = e~ ~ . 

As a simple check of our calculation, note that if one sends 7 to (keeping ( fixed), 
since the constraint (5.5), which was the only difference with the disordered phase, disap- 
pears, one should recover the results of the previous section. This is indeed what happens 
when one replaces the theta functions with their q — > limit. Also, (5.22) has been 
numerically checked with high accuracy. 
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This concludes the calculation of the bulk free energy in the anti-ferroelectric phase. 
Restated more explicitly, this is the result we have obtained: the partition function 
of the six- vertex model on a N x N lattice with DWBC and Boltzmann weights given by 
(5.1) has the following large N behavior: 

lim Z]l^ = sinh 7 (1 - C) sinh 7(1 + Q — -^2- (5.23) 
tv^oo N ,v sy ,v s; 27 6» 2 (^) 

where ( = i/7, and the elliptic nome of the theta functions is q = . Note that this 
expression is different from the corresponding expression for PBC. Let us now consider the 
two limits 7^0 and 7 — > 00. In both cases we shall assume that ( remains fixed. 

5.2. Small 7 limit 

As one sends 7 to 0, one reaches the line of the disordered/anti- ferroelectric phase 
transition. As noted earlier, the bulk free energy of the disordered phase is essentially 
obtained from that of the anti-ferroelectric phase by setting q = in the theta functions 

7T 2 

(and performing the rotation 7 — > ry, t — > it in the prefactors). Considering that q = e~ 2? , 
we expect a very smooth phase transition. More explicitly, we have the following expansion 
of/: 

+ 2 E ^r^ 1 - cos(m7rt/ 7 )) (5.24) 

m=l 

After subtraction of the analytic continuation of the disordered phase free energy (note that 
this analytic continuation is trivial since / only depends on t/7), we obtain the singular 
part of the free energy, which has a leading singularity 

/sing = 4e- 7r2 /7 cos 2^ + ... (525) 

This is the same type of singularity that appears in the model with periodic boundary 
conditions [5]. In more physical terms, if we introduce a temperature T which is near the 
critical temperature T c , we have 

/ si ng oc e ~ c / VT ^ T (5.26) 
that is an infinite order phase transition. 



/ = log 



TV 



27cos(f^ 
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5.3. Large 7 limit 

Next, let us consider the 7 — > 00 limit, i.e. A = — cosh(27) — > —00. This is a typical 
zero temperature limit, and we expect that the free energy will be dominated by the 
contribution of a ground state. After a modular transformation, the bulk free energy reads 



F = — log(sinh(7 - t) sinh( 7 + t))-f 



t 



-I -logsinh( 7 + t)+t + 2^ 



x I e -2m 7 



27 ^ m sinh(2m7) 

m=l v ' 

We can interpret the first terms when 7 — > 00 



sinh 2 (m(7 — t)) 



(5.27) 



(5.28) 



as coming from the family of ground states described by Fig. 3. The pattern of a rectangle 
inscribed inside a square is reminiscent of the circle inscribed inside a square characteristic 
of the disordered phase [2] . 



a 



b 



/ \ 



2 / 



1 s ' 
1 \ / 



a. 



Fig. 3: Ground states of the anti- ferroelectric phase. In regions a and b the 
arrows are aligned, whereas in region c they alternate in direction. 



5.4- Subdominant corrections 

As a final note, it is interesting to understand why the approach of [1] fails in the anti- 
ferroelectric phase. There, the idea was to find an appropriate Ansatz on the asymptotic 
behavior of the determinant tn and plug it in the Toda equation (2.5). The simplest 
assumption is that only the leading behavior (bulk free energy) must be taken into account, 
which leads to replacing tn with cn e N ^, where cat = (Yln=o n ^) 2 - ^ ne Toda equation 
then reduces to the ordinary differential equation for /: 

/" = e 2/ (5.29) 
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We can now use some insight from matrix models to understand whether this as- 
sumption was justified or not. In the ferroelectric and disordered phases, we reduced the 
computation of tn to a matrix model with eigenvalues in one single interval [a, b] (disre- 
garding the saturated region which plays no role here); it is known that such models have 
a regular large N limit. In fact, in the ferroelectric phase one can easily prove that 

T N ~c N e N2 fe N «-W (5.30) 

up to only exponentially small corrections; whereas in the disordered phase, one expects 
an asymptotic expansion which starts with 

r N ~ c N e N * f N R C (5.31) 

and continues with inverse powers of N (note that this is not quite the usual topological 
expansion of 2D gravity since the potential is not polynomial) . In either case, the assump- 
tion on the corrections is valid, and indeed, one can check that the expressions (3.15) and 
(4.12) do satisfy the ODE (5.29). 

On the contrary, in the anti-ferroelectric regime, we have found that the support of 
the eigenvalues contains two intervals [a, a'} and [&', b] and therefore we expect to be in 
a situation similar to what was studied in [19,20]. The analysis shows that tn should in 
this case display a pseudo-periodic behavior, which is indeed what is found in numerical 
computations. More precisely, after some calculations along the lines of [20] , one finds that 



tat ~ cn 



N 

* gUQ) 

27 2 (^) 



6 



(|(l + C)tf)c (5-32) 



■2 



where £ = t/7 and the elliptic nome q of the theta function is as before q = . The 
constant C depends only on 7. One can check that the right hand side of Eq. (5.32) does 
satisfy the Toda equation (2.5), even though the bulk free energy alone does not satisfy 
the ODE (5.29). 
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